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ABSTRACT 


In the paper, it is constructed with difference schemes which approximately mixed problem for heat equation and shown 
their stability. There exist various methods to develop difference schemes which are mainly based on exchanging 


derivatives with difference schemes. 
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INTRODUCTION 

Preliminaries 

Problem 1: We consider the following mixed problem for the equation below 



0 ) 


With boundary conditions for X = 0, (f, y) G R + 


u t ~a { u x -b l u y = 0, 


( 2 ) 


For y = 0, ( t , x) G Rf 


u t -a 2 u x —b 2 u y = 0, 


(3) 


And initial condition for t = 0, (x, y) G Rj 
u( 0, x, y) = <p(x, y), u, (0, a, y) = ip(x, y) 
In the domain Rl = { ( t, X, y)| t , X, y > 0}. 


(4) 


Here a, , h , , b 0 E R and r —■> 0 


1, , u-2 9 U 2 
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If for the conditions (2), (3) Shapiro-Lopatinski condition holds [17,50] 

• X = 0 in a, > 0, \b x | < 1; 

• y = 0 in a 2 > 0, \b 2 \ < 1; 

We rewrite the problem (1)—(4) in new coordinates system C,. 8 (x = r cos 8, y = r sin 8, £ = In r ) 
e 2 ^u tt -u ee —Ug =0, for t > 0, 0 < 6 < y, R l 


A 7* 

e^u, +a l u e -b t Uc = 0, for 6 - — , t > 0, ^ e R ] 


e^u t -a 2 u e -b 2 u s = 0. for 8-0, t>0, <%eR l 


(5) 

(6) 

(7) 


Then boundary conditions, t — 0, 0 < 8 < — 

m| = ) = (pie* cos 8,e^ sin#) 

u t | = y/(i;,8) = y/( e ^ cos 8,e^ sin#) 

In |^| —> oo an initial conditions u t = o(e ™), u g —o(e ■*”), u^—o(e Obtained problem (5)-(8) in the 

domain / > q o <6 < — fe A’ 1 is reduced to the following mixed problem consisting of system of symmetric f-hyperbolic 
’ 2 

equations 


f £ d d d ) 


r\-"‘»- c ^ +a r o 

(9) 

Boundary conditions 


# = y, f> 0, R' for +a x & 2 -bfy - 0 , 

(10) 

8 = 0, t >0, % e R 1 for - a 2 8 2 - b 2 8 2 = 0 , 

(11) 


Initial condition, for t = 0, 0 <8<—,^gR V — (c s lj/{8,^), (p g {8,£), (p^{8,^)} and when 
|£|—V = o(e^) ( 12 ) 

Here 
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according to [8], taking notation V = e 




Y 


, we get equality of dissipative energy integral 


et(A 0 Y,Y) t -(B 0 Y,Y) g -(C 0 Y,Y)t + ([Q 0 +Ql-2<uC 0 + ^-B 0 ]Y,Y) = 0, 

Clu 


(13) 


And assuming // = — aS We ^ aS w ^ en 1^1 °° IMI = ^~> 0, we integrate it in the domain 

n = j(#,£)| e R, o < & < 


d_ 

dt 


\\\e S (\Y,Y)d$dey\{(B„Y,Y\ 




rr 

d 

\ 

+ ff 

Q 0 + Qo ~C 0 + —— B 0 

Y,Y 

JJ 1 

77 V 

d9 

) 


d^d9 = 0 


a m=To 

B o (0) = T* 


We use the following equalities (see [3],): 

f H{ 6 ) O ^ 

v 

f O 
H(0) 


O H(9)) T °’ 
-H( 9 )' 


v 


c 


O 


T 

1 0 ’ 


y 


C o (0) = T o 


H(0) O 


\ 


O 


H(0) 


T ■ 

1 0 ’ 


Qo(0)+Q»(V) = 


2m(9) 0 k(0) 
0 0 0 
v k(0) 0 0 


= T 


~H{0) L(0) 
L(9) H(0) 


0 ’ 


O 


A) (A 

v 0 0 y 


, L( 9 ) 


-1(0) m(0) 

m(9) 1(9) 
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H(0) 


k(0)-m(0) -1(9) 

-1(0) k(9) + m(9) 


V2 


1 0 -1 

0-10 
0-10 
1 0 1 


Now we study ^(AqFjF) and 


Qtt + Qu -C 0 + B 0 
a 9 


Y Y ■ I* 1 order that this inequality 


holds e s (A 0 F, F) = {A 0 V, V ) = (HW, ,Wj+ ( HW 2 , IF,) > 0, it has to be H > 0, where 


W = 






w, 




V w 4y 


= T 0 V 


If k = k(0) > 0, k 2 (0) - m 1 (0) -1 2 (0) > 0 ,thentf >0. 


rr d 

\ 

f 

7 -H 

L ) 


f-H 

o' 


f 

O 

-H'' 

\ 

Qo+Qo~ C o + —B 0 

^ cl u 

Y,Y 

= e ~ 1 



— 



+ 




v,v 

/ 

V 

[ L 






V 

-H’ 

o 

J 


If H' =- H — L , that’s if (k — m)' = —l, (k + m)' = l, —l' = m, then square form is equal to 0. 

d9 

Finding the solution of these differential equations, we have 


k(0) = 


k( 0) 

f 

k 


n y m=< 


l (0) cos 9 - m(0) sin 0 


m| y |cos# + / 


n \ . „ , m(0) = 


— | sin 

v 2 


m( 0) cos 9 + 1(0) sin 9 


-1 


( Jl 

— \cos9+m 

v2 


f - | sin 9. 

v 2 


Now we 


consider square forms — (B q Y ,Y)\ g _n and (B q Y,Y )| . 


We rewrite Boundary 

(B 0 V ; V)\^ = e-t 


7t 


condition (10) as W t = SW 2 when 0 = —, where 


5 = 


2a t 1 -b t 
\+b t 
0 


1 + h 
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With the help of this equality, we get 

- {B 0 Y,Y\^ = e-^H + Hs]v 2 ,W 2 ] gJL 

And analogously for 0 — 0 rewriting boundary condition (11) as W l = RW-, we find 




0=0 


Where 

^ 2a, 


R = 


l-/z, 


1 + b, 1 + b. 


1 


0 


For n +1 solutions are obtained by a formula explicitly while the solutions for n of the difference schemes are 
known, this scheme is called explicit. Even though samples of some difference schemes to be considered look like the 
samples of explicit schemes, they are not actually. 


Explicit Right Difference Scheme: To solve problem-1 numerically, we employ explicit right difference scheme which 
approximates differential problem. For this, we rewrite the system (9) in the following form (7): 


w_a[Ml_ c w + 

^ dt dd °d£ 


Qo RC q + B 0 

a 6 


Y = 0, 


e% > f - - 5 o tr C ° l 7 + [e ° “ MC ° = 0 


(14) 


(15) 


We multiply the systems (14) and (15) by D — didg^y^, y^, y 3 ) from left side. Adding obtained systems, we find 


Ay 

CD\ — + e^D\ D - L r°; J - DB„ -J-- DC, 

dt 


37 n 3[B 0 T] 

dt 


de 


37 
J °3 e 


37 
0 3£ 


(16) 


- DC 0 + D[Q 0 - pC 0 + A B]Y + D[Q 0 - M C 0 ]7 = 0. 

3 £ dd 


In the domain? > 0, (O, TJ , we employ mesh with step-sizes respectively 


A? = A,, A<9 = A e , A£ = A^. 

We take the following notations and define norm as follows: 

Y" - Y(n A,,/ A e ,j A^) = (y^n A t ,i A g , j A f ), y 2 (n A t ,iA g , j A f ), y 3 (n A t ,i A g , j A,))', 
( = 0,/, n, 1/1 = 0,1,...; 



:A„A t y y 




2/ 

1=0 j=- 
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L - (1,1, = 


Using the above notations, we develop difference scheme approximating (16): 


^G 

y n _ v n Y n — Y n Y n — Y n 

- D? +lj ( B„), +1 Mi ij -DZ(C 0 l^±-^--D?JC 0 l + 

Ag Ag 


+d; 


2Qo 2//C 0 + B 0 

eld 


Y" = 0, i = 0,1-1, hi = 0,1,2,..,« = 0,IV-1, 


for Z = 0, | j| = 0,1,2,... 


(yi)S;-^(.v 2 )o; -b 2 (y 3 r 0J = 0, 

for i = I, \j\ = 0,1,2,... 


(17) 


(18) 


(y l r l j+a l (y 2 r Ij -b l (y 3 ): J = o. 


(19) 


for n = 0, z = 0,l,...,/, |yj =0,1,2,... 

/ 

^ = (e*V(£.,3), ^'#(£.,3), . 


( 20 ) 


Sample of this difference scheme as shown in figure-1.1, consists of the systems of equations which are not linear. 
It is easily seen that, this scheme approximates differential equation with first order. The boundary conditions are precisely 
approximated. As we mentioned above to compute schemes approximation error exact solution is calculated by the scheme 

(17) and we denote this error 8f h as the norm of the vector E ( t n . 0 j . ). Approximation of difference scheme for the 
sample equation is shown in the figures 2-4 and table-1. Here t w = 0.3, t 20 = 0.6, t m = 0.9, 6 1 = 0.524, = 3.2, 

max(,E I0 ) = 7.015x 1CT 7 > max (E 20 ) = 5.729xl(T 7 , max(£ 30 ) = 4.719x10 7 . 

i,j hj l ’J 



On’ ^i+l 

Figure 1: Difference Scheme. 
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Figure 2: Values At t 20 = 0.6, 0, = 0.524 of the 
Error (E n ). of the Sample Problem. 



Figure 3: Values At d 2 = 0.524,4=3.2. 



Figure 4:. Values At t 20 = 0.6,4 = 3.2 of the Error 
(E^ofthe Error (£) ;/ 


Table 1: Value At t 30 = 0.9 ( e „) °f Approximation Error 



1=0 

1=1 

1=2 

1=3 

1=4 

1=5 

j=0 

4.719e-7 

1.91 le-7 

6.868e-8 

1.762e-8 

1.244e-7 

3.42e-7 

j=l 

3.457e-7 

1.367e-7 

4.735e-8 

1.546e-8 

9.432e-8 

2.455e-7 

j=2 

2.482e-7 

9.721e-8 

3.410e-8 

9.628e-9 

6.457e-8 

1.734e-7 

j=3 

1.755e-7 

6.815e-8 

2.416e-8 

5.942e-9 

4.378e-8 

1.208e-7 

j=4 

1.224e-7 

4.718e-8 

1.687e-8 

3.638e-9 

2.943e-8 

8.310e-8 

j=5 

8.437e-8 

3.231e-8 

1.164e-8 

2.213e-9 

1.963e-8 

5.659e-8 

j=6 

5.754e-8 

2.191e-8 

7.942e-9 

1.337e-9 

1.301e-8 

3.818e-8 

j=7 

3.889e-8 

1.473e-8 

5.369e-9 

8.037e-10 

8.572e-9 

2.555e-8 

j=8 

2.607e-8 

9.832e-9 

3.601e-9 

4.802e-10 

5.616e-9 

1.698e-8 

j=9 

1.955e-8 

7.938e-9 

3.517e-9 

7.528e-10 

2.541e-9 

9.788e-9 
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Theorem 1: If for (18)—(19) it holds Shapiro-Lopatinski condition, that’s if « > 0, |fc,| < 1 any > 0, |fo,| < 1, then difference 
scheme (17)-(20) is stable with respect to sj J" energetic norm, here j„ _ A A ^ X ( A' / V)" ’ 


Proof. We multiply (17) by vector L and for convenience, we perform this operation for each variable separately: 


A (A), 


Y n+1 — V n A ^ 

-2-—,L 


y 




y«+i _ y" ^ 

-—,L 


v 




+ i ; ( ^ +I - ^), (^)" +1 )= i ; ([ ^ (^ r 1 - , [ a ,^]" +1 )= 

=^{((A,),(r 1 -^)^;1+((A ) ),(r 1 -^)^r)}= 

= j - {((A)), r .1®)- ((a ), y ; a; )+ ((a ), *r‘. *r‘)- ((a ), y ;. at 1 )}= 

= ^ ((a), y;r . Y t )- ^ ((a), y ;; . a ) = ^ ( aa (a a +);; 


Here we used D '-L = K”, D" 1 L = F.” +1 and A n = A n . 


7 n +1 


y y 7 y 


D" 


. IvL -Wi 


,L 


+ 


1 .L 


r(lBoA'„,-[Bor]",i:; i )+ 


^9 

- (M m y ;, y,: tJ ))=A(s. y , ry M1 - - Av. y% : 


V" _V' 

d;(c 0 \ g+1 A ,l 


+ 


D"{C 0 \^-^,L 


A, 


= A- (((Co), i£,, D -l)- ((c q ), r ;, d;l)+ 


+ «C 0 ), i£,. o,;„c)-((c„), r;, z.)) = +(c 0 r, r);„ - A (c„r, r)’ ; 




20o c o+ m A) 

dd 

2&-C„+A b o 

dd 


Y" I 

y ’ L 


7 V 
A 7 


y« y« 

1 ij U 


VL 


20,, “ A) + “ A) 

dd 

00+00 “ A) + ~ A, 

dd 


y^d^l 


yn yn 

y ’ x ij 


here it is employed fe 0 ^, i") = y ((-2 0 ^",^)+ fo#.1?))• 

Summing up above equalities, we get discreet analogue of differential representation of energy integral: 
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Go + Go ~ Gi + — B u 
cW 


Y 


( 21 ) 


Y,Y 


= 0 


Multiplying both sides of (21) by At, , we sum up in i from 0 to / — 1, in j from — to + °° and using 


Y n 


(wf —> 0 when \£\ —^ 00 , we have 


II II? II M? 

y" +1 _ y" M " 


"A> 


ll/io 


= A, ■ A { ■ £ {(s„r,r)''. -(Bo s ^ 

y'=-“ i=o ;= 

One can easily check that the equalities hold: 


G„+Go-c 0 +—5 0 
dO 


Y 


Y,Y 


am)=t; 


\ 


H(0,) O 
V o H(0,), 

-Hid,) O 


f 


T 0 > B o (0i) = T* 


O 


-Hid,) 


c o (0,)=t;\ ' K'o=\ °|> L (3) = 

0 ' 01 <3 (<9,)J 0 1,0 0 


-HO,) O 
r ~W) Y6,) 


T 

-'ll ’ 


m(e t ) 1(0,) 


' Hm J k W- m W ~ l W s 

-1(0,) k(0,) + m(0 


V 


6 0 (^) + 6o(^) = 


2 m(0,) 0 k(0,) 
0 0 0 
v *(3) 0 0 


A 


-T 

± r\ 


-H(0 ;.) HO.) 
v *■($) #(3) 


r 

■* 0 ’ 


T = _ 

° V2 


1 0 -1 

0-10 
0-10 
1 0 1 


0, =iA g ,i = 0,1,...,/ • 


According to the above 

^'((AAW) =‘’H\r.r): = <‘“'*'(4 X v); = 




if A:(f? ) >0,k (6, ) — m (6, ) — I “ (6 ,) > 0, i = 0, 1 , then // (# ; ) > 0 and hence 

/((a,)^;^;)>o, ( 22 ) 


where . 

If 
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d_ 

d0 

then 


r k'(0)-m(0) - l'(0) 

-l'(0) k'(0) + m(6) 


me) = 


-1(0) m(0) 

m(6) 1(0) 


= L(0) 


Q 0 +Qq-C 0 +—b 0 
dO 


\ 

Y"Y" 

/ 

'(-me/ 

my 

y ’ y 

-i J 

V 

\ W) 

me/, 


(23) 


f -H(G/ O ) f O -H/qy 


o HiBy 


-h\b, ;.) o 


w;,w;; 


= o. 


Hence we find (k(0)-m(0)) =-1(0), ( k(0) + tn(0 )) =1(0), -l\0) = m(0) 

k'(6) = 0, l\d) = -m(0), m(0) = 1(0) ■ Solving these differential equations, we obtain 


k(0,) = 


\k(0o) 

iw 


m= 


\l(0») cos 0, -m ( 0 O ) sin 0 j 
I in ( 0 I )cos 0, +1(0/) sin 0 i 


, 111 ( 0 / = 


f m ( 0 O ) cos 0 t +1 ( 0 Q ) sin 0 { 
\-l(0/) cos 0 j +m(0 l )mO i . 


-(b u y,yT Ij =-((B 0 ),e %;;.)= 


= —e 


f r O -H( f)) ^ 

W 0 V,V 


V°l 


=^ lm®w 2 ,wX.+{ h($w v w 2 X. } 


'ij 


~H( f) O j 

We rewrite boundary condition (19) as (W^ . = S(lF 2 )" ■, where 

r 2 a, 1-bA 


S = 


1 + b { 1 + b, 
1 0 


Due to this equality, it holds 

-(S 0 F,F);.=^([s* J f/(f) + //(f)5]W 2 ,W 2 )\ 


According to lemma D.L. Tkachev, M.V. Gomolina ([7]) if for boundary conditions Cl^ > 0 and |Z?j | < 1, 

[s - H(f)+ff(f)sl .,>0 


From this it yields 

~(B 0 Y,Y)"> 0. 


(24) 


Analogously we rewrite C b 0 y,yX taking into account that boundary condition (18) is (W,)”. =R(W 7 )" )j as 


o j 


'0j 


{B 0 Y,Y) n 0j =-e^ ([R*H(0) + H(0)R-\W 2 ,W 2 )'‘ oj . 
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where 


R = 


2 a 2 1 - b 2 

1 + b 2 1 +b 2 

1 0 


. If a-, > 0 and b-, < I for boundary conditions, then according to Lyuapunov 
theorem [/?*//(0) + //(0)/?] ;= o < 0 ■ Consequently, it yields 




(25) 


According to (22)-(25), we have the following energetic estimate 
2 

A) 


Y > ,+ 1 

2 

< 

y n 


A) 



From Y" = e^'V" it is easy to understand that 

v v 


a az z Awr= A AZ z ^ (ka i 4 j vr ^ iy r )= 


i=0 y=— 

/—I oo 


1=0 ;'= 


= a az z (k)a +i a ,+i )= a az z (A,v.v)r 

z'=0 j=—°° i= 0 j=—°° 

Therefore, it holds 


r = A AZ ZAV,A< A AZ Z(A,v,A=^° 

1=0 j=—°° i =o j——°° 

Theorem is proved. 

To solve the first problem, we make use the following difference scheme: 


u;r-u; 


A)"=o. 


(26) 


(27) 


The above difference scheme is obtained using U t ~U t = 0 . Now we multiply (27) by the vector 2 U"f : 


2 2 (UJ2 a , (((/,)".,(/") = 0 


and from this equality, we have 

(£/;*',£/“ l )-(t/;,£/;)-A 1 ((£/ t );,(£/,);)-A r (c/" 1 ,c/")so. 

Here Cauchy-Bunyakovsky 2 (U,V) <([/,[/) + (V,V) is employed. Obtained inequality is multiplied by 
A 0 A c and is summed up with respect i from 0 till / — 1, with respect to j from —°o till +°o : 
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>'=0 j =-°° i =0 j =-~ 

+^,4,4,21; 


i=0 ;=-■ 

We rewrite above inequality 


U' 




< \\U n \\ +A, 


(u t y 


A) 


+ A, \\U' 


llA) 


Writing (26) as 

c, {a A £ £ (u„ur + (u x ,u,C + (u„u,r\<r< 


=0 j=- 
I -1 


<C 2 A A x Z (.u„u l y i . + {u x ,u x Y ij +(u y ,u y ) n 

i =0 j=-°° 


we add it to (28). It yields discreet analogue of (0.10) 


u 


n +1 


< const 


U n 


(28) 


where 


U' 


= aaI z Wn 


i=0 >=- 
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